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We present a theoretical study ol the orbital magnetoresistance in a unixial anisotropic metal 
within the relaxation-time approximation. The appearance of a new dimensionless scale, 8 = 4t±/ep, 
allows the possibility of a new region at intermediate fields where the magnetoresistance is linear in 
applied magnetic field for currents flowing along the unixial direction. (Here, t± characterizes the 
bandwidth along the unixial direction.) In the limit of large anisotropy (small 5), corresponding to 
a quasi-two-dimensional metal made up of weakly coupled layers, we obtain an analytic expression 
for the magnetoresistance valid for all magnetic fields. We test our analytic results numerically and 
we compare our expressions with the c-axis magnetoresistance of S^RuCU. 



I. INTRODUCTION 

A growing number of compounds have been syn- 
thesized whose crystal structure consists of weakly 
coupled metallic layers. Foremost among these are 
the cuprate metals where the two-dimensional na- 
ture electronic structure has provoked much specula- 
tion about the nature of the resulting metallic state. 
The issue of whether in-plane excitations can move 
coherently between copper oxide planes remains con- 
tentiousa. Several other metallic compounds with a lay- 
ered structure — for example organic compounds based 
on the bis-ethylcnedithiotetrathiafulvalene molecule and 
other metallic oxides such as Sr 2 RuC"4 — do possess a 
Fermi surface which is shaped like a slightly warped 
cylindcrcl. In this paper we study magnetotransport in 
a quasi-two-dimensional (2D) metal in order to provide 
a benchmark against which more exotic types of behav- 
ior can be compared. Specifically we study out-of-plane 
transport using the relaxation-time approximation in the 
presence of an in-plane magnetic field. Our main result 
is that there is a "Kapitza" regional where the trans- 
verse magnetoresistance Ap c / p c is proportional to the 
applied magnetic field. We justify this with an analytic 
expression for the magnetoresistance in the quasi-two- 
dimensional limit which is valid for arbitrary magnetic 
fields. We also obtain an expression for the magnetore- 
sistance with any degree of unixial anisotropy that can 
be evaluated numerically. 

The study of magnetotransport within the Boltzmann 
formalism is long established and is well described in 
a number of classic textaJu. Interpreting magnetore- 
sistance measurements is complicated by the fact that 
the magnetoresistance is identically zero for an isotropic 
metal. The amount of anisotropy determines the mea- 
sured magnetoresistance, which is therefore rather sensi- 
tive to the detailed properties of the material. Analytic 
results are usually limited to the very weak or very high- 
magnetic-field regimes. At low fields the Zener- Jones 
expansion yields a magnetoresistance quadratic in the 



magnetic field with a coefficient depending on the vari- 
ation of the mean-free-path around the Fermi surface. 
At high fields the magnetoresistance saturates when cur- 
rent flow is along closed Fermi-surface directions or main- 
tains a quadratic field dependence for currents along open 
Fermi-surface directions (see page 118 of Ref. ||). To our 
knowledge, there have been no analytic expressions for 
the magnetoresistance of a realistic bandstructure which 
interpolate between these known limits. 

In this paper we present a calculation which, while re- 
specting the high- and low-field results mentioned above, 
also applies at intermediate fields where we find a linear 
magnetoresistance. We have obtained an analytic expres- 
sion for the magnetoresistance, valid in the limit of strong 
anisotropy which we believe to be the first straightfor- 
ward example of a magnetoresistance formula valid for all 
magnetic fields. This result should prove useful in char- 
acterizing the properties of quasi-two-dimensional metals 
using transport measurements. The paper is organized as 
follows. We first present a simple calculation of the mag- 
netoconductance that illustrates how having a warped 
cylindrical Fermi surface can give rise to a linear mag- 
netoresistance. We then present a more formal solution 
of the Boltzmann equation and derive the conductivity 
tensor for arbitrary levels of anisotropy. Our analytic re- 
sult emerges as a limiting case. Finally we discuss the 
implications of these results for experiment and compare 
with known data. 



II. SIMPLE PICTURE 

We consider a metal with the following dispersion re- 
lation 

h 2 

e(k) = — (kl + k 2 y ) - 2t± cos (k z c) + 2t± . (1) 

We have adopted the customary notation that directions 
in reciprocal space are labeled x, y, and z while the cor- 
responding real space lattice is defined by the a, b and 
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FIG. 1: The field geometry and approximate quasiparticle 
orbits — here shown for S — 2/3. (a) The magnetic field lies 
in the plane and the current flows out of the plane. The 
linear magnetoresistance regime occurs when the magnetic 
flux through a loop defined by the in-plane mean- free-path, 
la, and the plane separation, c, becomes of the order of a 
flux quantum or greater, (b) In the limit of large anisotropy, 
we can approximate the quasiparticle orbits as lines on the 
Fermi surface with constant radial angle 9 and constant in- 
plane velocity v a . This gives an approximate derivation of the 
low- field (B 2 ) and intermediate- field (|B|) regions. 



c directions. It describes free particles in the ab plane 
coupled by a perpendicular transfer integral, t± in the c 
direction to adjacent planes. The magnitude of c gives 
the spacing between planes that can be combined with 
t± to form an effective mass for out-of-plane motion 



mj_ = h 2 /2t±c 2 



(2) 



(This is the z-axis band mass in the limit < fx when 
the Fermi surface forms a closed spheroid.) 

We begin with an approximate derivation of the new 
regime by considering the magnetoconductance to lowest 
order in t±. Chambers' expression^ for the components 
of the conductivity tensor in a magnetic field within the 
relaxation-time approximation is 
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dS 
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VityvAtye-^dt. (3) 



For each area element of the Fermi surface, dS, we in- 
tegrate the velocity v(t) measured along a semiclassi- 
cal quasiparticle orbit. These orbits are defined by the 
Lorentz equation of motion 



n — = — ev x B , 
dt 



(4) 



where v = Vfce(fc)/7i. In this paper we will be interested 
in configurations where the magnetic field is parallel to 
the planes and the current is flowing perpendicular to the 
planes [see Fig. (Ha)]. 



To lowest order in t± we can ignore both the z-axis dis- 
persion in the equation of motion and the closed orbits. 
The rate of change of k then only depends on the an- 
gle 9 between the magnetic field and the in-plane Fermi 
velocity of an electron [see Fig. (0b)] . So we have 



. dk z 

l — ; — 

dt 



-ev n B sin 9 



(5) 



where the in-plane Fermi velocity, v , is constant here. 
Hence we see that 



ck z (t) = ck z (0) + tt c ti 



(6) 



We can determine how the velocity in the c direction 
changes as the quasiparticle moves across the Fermi sur- 
face by combining the above result with Eq. (uh giving 



v c (t) 



2t±c 



sin [k z c + Q c t sin 9] 



Here we have defined a "cyclotron" frequency 



fir 



ev a Bc 



(7) 



(8) 



which is the fastest rate at which quasiparticles traverse 
the Brillouin zone. It turns out that this sets the scale of 
the crossover from weak (fi c T <C 1) to intermediate fields 

(fl C T>l). 

Substituting this velocity, v c (t), into the expression for 
the conductivity Eq. (|J), we may integrate over time and 
k z to give the out-of-plane conductivity expressed as an 
integral over the orbits, namely, 
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(9) 



Integrating this is straightforward and yields the follow- 
ing magnetoconductance 
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Aa c {B) = 

<J C {0) y/T+Sfy 



(10) 



At low fields (O c r <§C 1) this gives the usual quadratic 
field dependence as obtained by the Zener- Jones expan- 
sion. However in the intermediate field regime (f2 c r > 1) 
the magnetoconductivity falls off as 1/B. This is the es- 
sential result of this paper. It arises because there is a 
range of cyclotron frequencies for traversing the Brillouin 
zone. Note too that the magnetoconductance becomes 
universal, independent of the degree of anisotropy and 
depending only on the in-plane properties. We can em- 
phasize this by writing the cross-over condition in terms 
of the magnetic length l 2 B — h/eB, namely fl c r = 1 cor- 
responds to 



w = 1 - 



(11) 



So the linear region is reached when approximately one 
flux quantum threads an area formed by the in-plane 



mean-free-path (l — v t) and the interplanar spacing. 
For a typical layered oxide (c ~ 12 A) one can expect to 
see this regime at 10 Tesla when the in-plane mean free 
path reaches around 500 A. As discussed later, experi- 
mentsQ on S^RuCU provide evidence for the validity of 
this expression. 

However, this cannot be the complete story since very 
general arguments show that a c must go as 1/ B 2 at high 
fieldso. Since the conductivity is the sum of the conduc- 
tivities from all orbits, the high-field form cannot be re- 
covered simply by including the contribution of the closed 
orbits: a 1/B 2 contribution from closed orbits will never 
dominate the 1/B from open orbits. The correct high 
field result emerges when we consider higher-order effects 
in t± for both the open and closed orbits in the exact so- 
lution of section III. We will see that the high field 1/B 2 
regime occurs when Q c t ^> yjep/t_\_. 

A further reason for a more detailed treatment is that 
the magnetoresistance is identically zero for an ellipsoidal 
Fermi surface with a constant r. This is because of a 
cancellation between the Hall and magnetoconductance. 
We therefore would like to verify that there is no such 
cancellation here. We will do this through an exact so- 
lution of the Boltzmann equation in which we compute 
all components of the conductivity tensor and hence the 
magnetoresistance. 



III. SOLVING THE BOLTZMANN EQUATION 

Our treatment will follow closely that of Abrikosovi. 
In the previous section we only treated the quasiparticle 
orbits approximately. Since the Lorentz force in Eq. (jij) 
acts perpendicular to the electron motion, energy is con- 
served and the electron is constrained to move along a 
constant energy line with fixed momentum in the B di- 
rection. Here we take B to be parallel to a. For finite 
t x , is no longer constant along the orbits and there are 
some closed orbits as illustrated in Fig. (0a). 

Quasiparticles at the Fermi surface determine the 
transport properties and we may identify two regions de- 
pending on the degree of anisotropy. We introduce a 
parameter which is a measure of this anisotropy 



(12) 



If 5 > 1, the Fermi surface is closed and the trajectory of 
all quasiparticles in momentum (and real) space follows 
closed loops. While the system may have an anisotropic 
effective mass, the qualitative features of transport will 
not be much modified from a typical three-dimensional 
metal. If S < 1 then we still have some closed orbits 
but there are now some trajectories that extend across 
the Brillouin zone in the c direction [see Fig. (§&)]. Only 
these orbits have been treated in section II (and then 
only approximately). We now aim for a more complete 
analysis. 
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FIG. 2: The exact orbits over the Fermi surface with 5 = 
2/3. (a) The quasiparticle moves along lines defined by the 
intersection of the Fermi surface and planes perpendicular 
to the field direction. These orbits can be either open or 
closed and are described by Jacobi elliptic functions, (b) We 
can label these orbits by a parameter a, the parameter of 
the Jacobi elliptic function. Here we illustrate orbits where 
a = 0.7,0.9,0.99,1.01,1.5 and 4. For 6 < a < 1 the orbits 
are open. For 1 < a < oo the orbits are closed. 



To consider the conductivity tensor we must solve the 
Boltzmann equation. Rather than use the momentum 
components k x , k y and k z , in the presence of a magnetic 
field it is more convenient to work in terms of a new 
coordinate system k x , e, and t. Here t is the time taken to 
move along the momentum orbits defined by the equation 
of motion, Eq. (|J). The advantage of this coordinate 
system is that the magnetic field is included implicitly 
and does not appear in the Boltzmann equation. Within 
the relaxation-time approximation, one may write the 
Boltzmann equation aaj, 

, * ,= . (13) 



Of 



= eE ■ v 



where the electron distribution function has been writ- 
ten as / = — ipdef 1 -^. This first-order differential 
equation may be solved straightforwardly. 
Substituting into Eq. (^) gives 

/;/,./; dk z 



dt 



eB 
\/a C 

UJ r 2 



2 m o 



fi 2 k 2 + Am t± [cos (k z c) — 1] 



dk 7 



where 



;sin 2 (k z c/2) 



8t±m 



— 2 9 " 

2m £F — h k 2 
We have defined a new "cyclotron" frequencyl 

eB 



eBc 2t ± 
h V m 



= n c \ - 

y/m mj_ V 4 



(15) 



(16) 



(17) 



4 



the fastest rate at which quasiparticles perform closed 
orbits. This is the natural scale for cyclotron motion 
perpendicular to the plane. As might be anticipated 
from a Bohr quantization picture, it also sets the scale for 
Landau-level quantization and hence the quantum effects 
which signal the breakdown of quasi-classical transport 
theory. The variable a labels each cyclotron orbit. It is 
bounded from below by 8 and we will use it to substitute 
for p a 



\/2m e — 8t±m /a . 



(18) 



The orbits are open for 5 < a < 1 and closed for 1 < a < 
oo. 

It is very unusual that Eq. (|l5|) is both integrable and 
its solution is invertible so that a closed form expression 
for the orbits may be found. The solution can be written 
in terms of the Jacobian elliptic functions 



k z (t) = 
v c (t) = 

v b (t) = 



— JacobiAmplitude 

c 



UJ r t 



, a 



2h 



JacobiCN 



LO c t 




LU c t 


—= : a 


JacobiSN 











2/7 



Cy/am±m 



-JacobiDN 



LO r t 



'a 



(19) 



(20) 
(21) 



These equations exactly describe the quasiparticle's mo- 
tion over the Fermi surface defined by Eq. (|l|) in the pres- 
ence of a magnetic field along the a direction. We have 
adopted the-aiotation of Mathematical] and Abromowitz 
and Stegun£2l in using the parameter a rather than the 
modulus, k — y/a, to define these functions. 

These periodic functions play the role of the trigono- 
metric functions that appear in the solution of the spher- 
ical problem. To make this more explicit, we can map the 
elliptic functions when the parameter a is greater than 
one (describing closedr-arbits) to those with a parameter 
(3 = 1/a less than onetl We may then write 

v c (t) = -^^JacobiDN[w c 7i,/3]JacobiSN[w c i,/3](22) 
cmj_ 



nm smlw c t) , 

(3^o c m± 



(23) 



2h 

Vb{t) = — < 



c y m±m^ 



- JacobiCN[w c t, 0] , 



,. 2h 

hm — 

/j^o c V mj_m 



■ cos(tu c t) , 



(24) 
(25) 



The limiting case of (3 — > describes all of the orbits 
when the Fermi surface becomes spheroidal (5 — > oo). 

To compute the conductivities we use the solution of 
the Boltzmann equation [Eq. (p"3])] and compute the cur- 
rent. This gives the Chambers' formula [Eq. (|^)] which 
may be written as (see Ref. j|) 



2e 2 B 
{2irhf 



dp a 



T(Pa) 

dt 



dt'vai^v^E^'-^ 



(26) 

This triple integral can be drastically simplified when we 
recall that the orbits are all periodic and so have a well- 
defined Fourier series. The FourierpSeries for the Jacobian 
elliptic functions are all tabulated!^. This allows one to 
do the integral over i! and, using the orthogonality of 
the components of the Fourier series, one can also do 
the integral over t. The algebra is somewhat tedious but 
the result is that the conductivity can be expressed as a 
rapidly convergent sum followed by a single integral over 

Pa- 

We will use the definition of the tensor conductivity 
ja = CapEp with the magnetic field along the a direc- 
tion. The only nonzero components of the conductivity 
tensor will be a aa , &bb, &cc and a oc = —a c b- There is 
no longitudinal magnetoresistance for the dispersion of 
Eq. (|l|) so a aa is unaffected by the magnetic field and, 
by rotational symmetry in the ab plane, will be equal to 
abb in the absence of a magnetic field. We can simplify 
some of the expressions by introducing a number of pa- 
rameters: a universal conductivity and the effective-mass 
ratio 



ctq 



2e 2 r 
c 3 mj_ 



(27) 



For the exact solution of the conductivity tensor we find the following components 



Obb 



min(l,(5) 



2irr<jn\f8 



mm(l,<5) 
1 



-da 



1(1,0 



1 + 8Si(a,u) c T/y/a) 
K(a)^Ja — 8 

5 2 (a,o; c r/- v /a) 



i(l,5) if(a) 3 a 3 V a — 8 



da + ; 



.da 



mm 



K(j3)y/1-P8dp , 



min( M -i) S 2 (P,L0 C T) 



n 



d/3 
d/3 



(28) 
(29) 
(30) 
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The remaining Hall component is given simply by 

eBr 

&bc 



m 



(31) 



In these expressions, K(a) is the elliptic integral and we 
have defined the following sums that involve, q, the nome 
of the elliptic integral q(a) = exp[— nK (1 — a)/K (a)] 



Si(a,x) 
S 2 (a,x) 



E 

n 

E 



q(a) n 



1 2 



1 + q(a) 2n 
nq(a) 



1 



1 + 



( mtx \ 
K(a)J 



1 2 



1 + q{a) 2n 



1 + 



( mrx \ 
\K(a)J 



, (32) 



(33) 



The summations above are over the positive integers (n — 
1, 2, • • • oo). The sums Si and S2 are the same as S\ and 
S2 respectively except they are summed over the positive 
half-integers (n = 1/2, 3/2, 5/2, ■ ■ ■ , 00). 

Equations ( f28|) to ( |33| ) are the exact solution for electri- 
cal transport within the relaxation-time approximation 
and are valid for arbitrary values of cf and tj_. To make 
further progress and to make contact with the result of 
the simple calculation outlined previously, we need to 
work in the limit of large anisotropy: <5 <SC 1. In this limit, 
the dominant term in the conductivity tensor comes from 
a cc [Eq. j30|)] and, in particular the small a range of the 
integral. It is therefore a good approximation to replace 
the integrand with its small a limit [i.e. take only the 
first term in the sum, and let K (a) — > tt/2]. Doing this 
gives the following approximate form for the conductance 



2tt 2 



arctan 



1 



(34) 

This result recovers the form of our approximate deriva- 
tion of the result [Eq. (|ic|)1 but remains correct in the 
extreme high field limit. With l/\/S 3> £I c t » 1 wc 
can replace arctan by tt/2 and we have the linear magne- 
toresistance regime as before. However, the correct l/£> 2 
asymptote is recovered when Q c t 3> X/VS. 

This is a better approximation than our previous treat- 
ment because we are giving an exact treatment of the 
lowest Fourier component of the open orbits. We are di- 
viding the conductivity into a sum over quasiparticle or- 
bits and now each of these becomes a Fourier series. The 
additive nature of the conductivity and the Fourier series 
means that each component must contain the physics of 
the high-field asymptotic limit as well as the linear inter- 
mediate field regime. We can compare the role of higher- 
order Fourier components which are neglected in deriv- 
ing Eq. ( p4[) by comparing with numerical treatment of 
equations |28| to [33| The dashed line in Fig. || shows the 
magnetoresistance keeping only the lowest order compo- 
nent [Eq. (34)] for S = 1CP 4 . The deviation from the 
numerically exact result indicates where high order com- 
ponents become important. However since all of these 
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FIG. 3: A numerical treatment of the magnetoresistance for 
arbitrary fields and Fermi surface topologies. These are com- 
puted using the results of Eqs. (J2s|) to (|3^). In the limit of 
S ^> 1 the magnetoresistance vanishes as expected in limit of 
a spheroidal Fermi surface. The previously known asymptotic 
regimes are unshaded [always B 2 at weak fields with a high- 
field asymptote that is B 2 for open Fermi surfaces (S < 1) 
or saturates for a closed Fermi surface (S > 1).] The lin- 
ear regime is heavily shaded while the crossovers are lightly 
shaded. The dashed line shows the result (for 5 = 10~ 4 ) from 
the lowest Fourier component of the open orbits [Eq. (BJ)]. 
Adding a numerical factor based on the asymptotic limit from 
all orbits [Eq. ( |35[ )] is also shown for this S but is indistin- 
guishable from the numerical result. Note that S < 10~ 2 for 
all three bands of S^RuCU. 



higher order terms contribute to the B 2 asymptote, they 
can be taken into account by including an extra numeri- 
cal factor in Eq. (34) so that the numerical and analytic 
results match in the high field, small 6 limit. This gives 
the following interpolating expression 



Ape 

Pc 



2arctan(l/^ + 0.0263 S{Q c t) 2 ) 



- 1 



(35) 



where the factor of 0.0263 is obtained by summing all or- 
bits numerically in the high field and small S limit. This 
function is plotted in Fig. ^ but is virtually indistinguish- 
able from the numerical result. The appearance of a new 
region where the magnetoresistance is linear in B may 
clearly be seen (Fig. ||) as the dispersion becomes more 
two dimensional. 

For completeness we consider also the Hall resistivity 
with current along c and the Hall voltage being developed 
along the b direction. We find 



PH 



'be 



eBr 



TOo (766 



(36) 
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2 4 6 8 10 12 14 

Magnetic Field (B) /Tesla 

FIG. 4: The transverse c-axis magnetoresistance of SnRuCU 
(data after Hussey et al, Ref. Q). This clearly shows the linear 
regime we have predicted. The development of the regime at 
relatively low values of Ap c /p c ~ 0.2 as opposed to ~ 1 is 
an indication that more than one band with differing mean- 
free-paths are involved in the c-axis conductivity. This is 
demonstrated by the fit, which includes a band with low c- 
axis conductivity and short mean-free-path. This is consistent 
with the band structure (Ref. ^) and the Hall coefficient 
(Ref. O). 



Thus there is no anomalous regime in the Hall resistivity, 
which remains linear at all magnetic fields and reflects the 
carrier concentration in the usual way. 



IV. COMPARISON WITH EXPERIMENT 

Sr2RuC>4 is probably the best characterized two- 
dimensional metal. The current experimental c-axis mag- 
netoresistanceu clearly shows the linear magnetoresis- 
tance regime we have discussed [see Fig. (^)j. We now 
discuss what quantitative information we may determine 
from this. 

Detailed de Haas van Alphen studies give a very clear 
picturp_of the degree of warping of the Fermi-surface 
sheetslij. These are expressed as variations in the ra- 
dius of the Fermi surface in the xy plane as the Brillouin 
zone is traversed along the z direction. To translate to 
our notation we note that 



Ak r 



2t x 
e F 



(37) 



In general, terms in the c-axis dispersion can involve 
cos(^fc z c) where v is an integer. Thus far we have con- 
sidered v = 1. Terms with v > 1 increase the c-axis 
conductivity which depends on the square of the c-axis 
velocity (~ v 2 ). Furthermore, there can also be an an- 
gular variation of the c-axis dispersion, tx, within the 
plane. This not only modifies the numerical prefactor in 
the conductivity but generally means that the c-axis mag- 
netoresistance becomes dependent on the orientation of 



the field within the ab plane. This is .known to be impor- 
tant in the cuprate superconductors!-! Finally, for the 
particular case of Sr 2 RuC>4, there are three bands which 
give additive contributions to the total conductivity. 

Because of the three bands and angular variation of 
t_i_, it is misleading to use the formulas we have de- 
veloped to determine a quantitative measure of Fermi- 
surface anisotropy from the magnetoresistance. Indeed 
we have argued that the magnetoresistance becomes uni- 
versal for each band in the 6^0 limit. Instead we can 
use the experiments as a probe of in-plane properties — 
the mean-free-path. 

To obtain a good fit to data while introducing a mini- 
mum of free parameters, we consider 

a cc (B) =0i+ ■= t - _ . (38) 



a\B 2 



This represents one fluid of electrons (f3\) with a short 
mean- free-path, I <C 1%/c, and therefore insensitive to 
the fields, and second (/3 2 ) with a much longer mean- 
free-path [following Eq. (MY)}. The magnetoresistance 
depends only on fii/fa ~ 0-1 an d oi\ ~ V0.04. This, 
is consistent with de Haas van Alphen measurementsEj 
which suggest that one Fermi surface sheet, 7, is con- 
siderably less dispersive in the c-direction than the other 
two. In addition, the assumption of a small mean-free- 
path on that sheet is also consistent with the Hall coef- 
ficient that remains strongly_temperature dependent in 
the regime of this experiments. a% may be related to the 
mean-free-path using l Q = aiU/ec which gives a value of 
l a ~ 2000A at 3K on the two sheets with the most z-axis 
dispersion. This is consistent with the observation of un- 
corjjcentional superconductivity in this sample at around 



V. CONCLUSION 

So to summarize: we have given an exact solution for 
electrical transport within a quasi 2D band structure. In 
doing so we find that the new dimensionless parameter 
8 = 4:t±/ep is important. This leads to a new region 
(V5 < uj c t < 1) in the magnetoresistance where the c- 
axis transverse magnetoresistance is large (Ap c /p c > 1) 
and linear in the applied field. An asymptotically exact 
expression for the magnetoresistance has been obtained 
in the limit of small S, i.e. the limit of weakly coupled 
2D planes. This new region has been observed at low 
temperatures in the quasi 2D metal Sr2Ru04. For the 
over-doped thallium cuprate there are signs that one is 
beyond the low field regime at 11 Tesla£3. 
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